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Massless QED in 2+1 dimensions (and in general for space-time dimensions 2 < d < 4) 
has very interesting features. It is not just super-renormalizable, it is ultraviolet finite. The 
usual perturbation expansion in the fine structure constant (which is now a dimensionful 
parameter) has severe infrared (iR) divergences which become worse as the number of loops 
increases. So QED3 provides an ideal platform for tackling iR divergences. This has led 
to an extensive study of this model pQ. It is clear that the usual perturbation theory can 
make sense only by some kind of resummation. In a 1/N expansion, N being the number of 
fermion flavours, there is a resummation of chains of one-loop vacuum polarization diagrams 
on every photon propagator. This changes the iR behaviour of the photon propagator from 
being inversely quadratic to inversely linear in momentum. Thus the iR divergence is 
softened. Even after this there are logarithmic iR divergences. (Throughout this paper we 
are concerned with the iR divergences in Green functions for non-exceptional Euclidean 
momenta and not the additional iR divergences for real processes.) The problem is to 
sum them up and extract the iR behaviour of the Green functions. We could handle this 
problem in the following way. We have shown in Ref. [2] that for a particular value (chosen 
to each order in 1/N) of the gauge parameter in a specific non-local gauge, the logarithmic 
iR divergences are absent. As a consequence the iR behaviour of the Green functions to all 
orders is known. The limiting behaviour is a conformal field theory where the photon has 
non-canonical scaling dimension one for the entire range of d, in contrast to the engineering 
dimension (d — 2)/2. The fermion continues to have the canonical dimension (d — l)/2. 

This behaviour for the fermions is of course gauge-dependent and special to this gauge. 
The Green functions involving only the photons is gauge invariant and the scaling dimension 
one would be valid in any gauge. Our specific choice of gauge has the advantage of extracting 
this information without being cluttered by the powers of logarithms in the intermediate 
stages of the calculations. 

Thus we know the iR behaviour of the Green functions to all orders for a particular 
choice of the gauge parameter in a specific non-local gauge. It is instructive to know how 
the logs add up for other values of gauge parameter and also in other gauges. In this paper, 
we predict the behaviour to all orders without detailed calculations. 

This is done using the Zumino identities [3J which exactly relate the Green functions 
in one choice of gauge to those in another. Such a relation is also called the LKF trans- 
formation this name has been mostly used for relation between various conventional 
covariant gauges (for example, between the Landau gauge and the Feynman gauge). We 
are interested in a more general class of gauges, including non-local choices. Although the 
relevant relation is contained in Zumino's paper we rederive it in a way which is suitable 
for our purpose. 

This far, our discussion has been restricted to a non-local gauge. It is also of interest to 
know how the Green functions behave in a conventional gauge such as the (local) Lorentz 
gauge. With our choice of the non-local gauge the longitudinal part of the propagator has 
the iR behaviour q^qv/q 3 , which is softer than that in the Lorentz gauge, q^,q v /q A - This 
feature led to only log iR divergences in fermionic Green functions which further could be 
cancelled to each order by adjusting the gauge parameter at that order. In contrast, in the 
Lorentz gauge, the iR divergences become increasingly worse with the number of loops. A 
resummation seems to be beyond reach. We find out the iR behaviour in this case also 
using the Zumino identities 
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Consider the generating functional of the Green functions in a general non-local gauge 
parametrized by function g(x,y) 



J exp (-\dA -g-dA + j-A + ri-ip + ri-ip) 

J exp {-\dA ■ g ■ DA) ' 1 J 

Here j stands for the measure 

J Vi)Vi)VA exp (-S) , (2) 

where S is the gauge-invariant Euclidean action. Also, is the source for the vector 
potential A^, and the Grassmann variables rj and fj are the sources for the fermions. We 
have used the notation 

OA = d»A\ (3) 

while the ■ operation signifies the inner product in the Hilbert space, involving integration 
over spacetime variables and summation of discrete (spin and/or internal) labels. For 
example, 



j . A = j d 3 x f{x)A»(x) , (4) 
dA ■ g ■ A = j d 3 xd 3 yd"A^x)g(x, y)d u A u {y) . (5) 
Varying the gauge function g by an infinitesimal amount 5g, we get 

5Z = -1 ((dA ■ 6g ■ dA) jtV>fl -Z(dA-5g- dA)) , (6) 

where 

/exp (-IdA-g-dA + j-A + fj-i/j + yfyO 
{ '1™ Jexp(-ldA-g-dA) {) 

is the expectation value of operator O in presence of sources j, r\ and fj. We take (O) 
without any suffix to denote the corresponding expectation value when the sources are 
set to zero. As a consequence of Ward identities, the correlations involving longitudinal 
photons in eq. © can be related to pure fermion correlations. This leads to the Zumino 
identities which we now derive in a form suitable for our purpose. 

Invariance of Z under the following change of integration variables (corresponding to 
an infinitesimal gauge transformation) in the numerator of eq. (0) 

8A»(x) = -d»e{x), , . 

6ip(x) = iee(x)ijj(x) , 8ip(x) = —iee(x)ip(x) 

gives the basic Ward identity 

{{d 2 g ■ dA + dj + ie (# - = (9) 
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where we have used the notation O x to mean 0(x). Also, d 2 g can be regarded as the 
product of the symmetric matrices (d 2 ) xy = S^(x — y)d 2 and g xy (note that eq. (JTJ) picks 
out the symmetric part of g). Applying the operation d^(8/ 'Sj' x (y)) on eq. (JUJ), we get 

(dAy (d 2 g ■ dA + dj + ic (# - . ^ = (d 2 J xy ) Z (10) 

where 5 xy stands for the Dirac-delta function 5^(x — y). Setting the sources to zero in the 
above equation, we get 

{dA y (d 2 g-dA) x ) = d% y . (11) 

This is the conventional Ward identity for the longitudinal part of the photon propagator. 
From the difference of eq. (jll)[) and (Z times) eq. ([lip, we get 

dAy (dA + {d 2 g y l ■ (dj + le (# -#)))) - Z (dA y dA x ) = 0. (12) 

In the last step, we multiplied with the appropriate element of the matrix (d 2 g)~ 1 , which 
is to be regarded as the inverse of the matrix d 2 g. Using eq. (fT2"j) in eq. (0), we get 

5Z= l -(dA-5g-{d 2 g y l -{dj + ie^-ri)))^. (13) 



3,V,V 



This is not yet the convenient form for our use. We proceed to eliminate dA from eq. 
(|13|) and arrive at an equation involving fermion correlations only. We apply the operation 
fj a (y)(S/6fj a (y)) on eq. © to get 

((&g-dA + dj + ie(ftf-7rf)) x (ftf) y } = -ieiifj^Jj^Sxy (14) 

and the operation r) a (y)(5 /Sr) a (y)) to get 

U&g-dA + dj + ie(fji>-ri$)) x (ri$) y ) _= ie ( S xy . (15) 

The difference of eq. f)14p and eq. (|T5j) yields 

^ (d 2 g ■ dA + dj + ie (fji[) - r]^)) x (fjif) - v4>)y). = ~ ie ( (# + ^) x ) j>v ^ <V (16) 

Multiplying with (d 2 g)~ x l and integrating over x, we obtain, 

- (dAg (fjif> - rtf) ) = (((d 2 gy 1 ■ (dj +ie(# -#))) (#-#)„) 

\ y / j',?7,77 \ \ / 2 "/ j,ri,r] 

+ie(d 2 g)- zy 1 ((fjiP+VTp) y ) jrj _- (17) 

Also from eq. 0, 

" (^>^fl = (( • + ie (#-#))) ), 
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Using eq. (fTSjl and eq. (|T7|) in eq. (|T3j) . and using the fact that (d 2 g) 1 is a symmetric 
matrix, we have finally 



52 = 2 



1 (dj + ie [jqip — ?70)) • 5 (d 2 gd 2 ) 1 • (dj + ie (t]-^ — 77^)) 



-e 2 6(d 2 gd 2 ) m ]( r] .^ + r ] ^)^_\. (19) 

Here we have used (d 2 g) 1 -5g-(d 2 g) 1 = — 5 (d 2 gd 2 ) , which can be easily checked by going 
over to the momentum space: (— k 2 g(k))~ 1 5g(k)(— kP'g^k))^ 1 = — {l/k 2 )5{l/g{k)){l/k 2 ). 
(We will come across explicit examples in momentum space later in this paper.) Also, 
(d 2 gd 2 ) Q Q appears in eq. (fTTfy as follows: 



/ 



d 3 x ((d'g)- 1 ■ 5g ■ (d 2 g) l )^{(^ + #) J 
d 3 x5 (d 2 gd 2 )~l((v^ + #)J 



-5{d 2 gd 2 )-l{{r 1 ^ + r 1 ^)) (20) 



since (<9 2 g<9 2 ) 1 depends only on the difference x — y for a translation-invariant gauge 



'xy 



function g. 

The dependence of all Green functions on the function g xy is contained in eq. Mty . The 
simplest case of the photon propagator is obtained by applying the operation S 2 / (§j^{x)b~j v (y)) 
on eq. (fT9"j) and then setting the sources to zero: 

5A^x,y)=d»d;6(d 2 gd 2 )- y 1 . (21) 

This is consistent with the Ward identity eq. fTTf . 

We now obtain the dependence of the fermion propagator on the choice of the gauge 
function g. Applying S 2 / (Sfj 1 (x)Sr]s(y)) to eq. (fT9"j) and setting all the sources to zero, we 
get 

8S r s(x,y) = -6F xy S lS (x,y) (22) 

F xy = e 2 {(d 2 gd 2 )^ - (d 2 gd 2 )-l) . (23) 

Integrating this equation, we relate the fermion propagator evaluated with two different 
gauge functions: 

S lS (x,y) = exjp[-(F-F°) xy \s° s (x,y). (24) 

Here S and S° stand for the fermion propagators in the gauges g and go respectively; F° is 
related to g° by eq. (J2H1). 

We have shown in Ref. |2| that if we choose a particular non-local gauge 



where F xy stands for 



(with p = Ne 2 /8) then it is possible to choose the gauge parameter a to each order in 
1/N such that there are no logarithmic corrections to the fermion propagator and other 
Green functions. As a consequence, for this particular value of gauge parameter (call it ao) 
the iR behaviour of the fermion propagator is that of the free theory with no anomalous 
dimension: 

S°»(x,y)~j^L. (26) 
\x — y\ 6 

The iR behaviour for other values of gauge parameter a within the same non-local gauge 
then follows from eq. (J24"j) and eq. (|25|). We now have F xy = af(x — y) where / is formally 
the matrix f xy = e 2 ((<9 2 (l + p/V^d 2 )d 2 ) 00 1 - (<9 2 (1 + p/y^d 2 )d 2 )~^). It is convenient to 
represent / by the Fourier transform 

r j3jL (l _ e ik.(x- y)\ 

«^»- 2 /^W' (27) 

[This is obtained by inserting the completeness relation for the momentum eigenstates in 
< x\(d 2 (l + yu/V— d 2 )d 2 )~ l \y >. Note that the factor of (27r) 3 in eq. (|27j) is consistent with 
that in the Fourier transform of \ xy = 5(x — y).\ We now get the propagator for the gauge 
parameter a: 

Ot — OtQ 



S " M =teJ Sa ' (X ' y) (28) 

where 

X(x,y)=exp[f(x-y)]. (29) 

The integral in eq. (|27|) is finite at both the ends k — > oo and k — > 0. For k — > 0, finiteness 
follows from fik ^> k 2 and exp(ik -(x — y)) ~ 1 + ik ■ (x — y) (actually, by symmetry, it is the 
0(k 2 ) term in the exponential which contributes). Now, as \x — y\ — > oo, k ■ (x — y) is no 
longer small, and the integral develops a logarithmic divergence as k — > 0. Thus, l/\x — y\ 
serves as an infrared cutoff for k, and for \x — y\ — > oo, we expect the integral to behave 
as Khx\x — y\ where k is a constant. This can be explicitly verified, and the constant of 
proportionality extracted, as follows. We have 

e 2 f°° dk / sin(k\x\) 



2tt 2 Jq k 2 + fj,k \ k\x\ 
e 2 f°° du ( , sinw 



27r 2 /i J pu 2 + u \ u 



(30) 



where p = l/(p\x\). Then, 



r df 
k = — hm p— 

p^o dp 



e 2 f°° dv ( sm(v/p) 

hm — — / -r 1 — p 

P-o 2vr 2 /i J (v + l) 2 V v 



-,2 



dv 4 



2ir 2 p J Q (v + l) 2 n 2 N 



(31) 
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which is finite and non-zero. This gives the iR (\x\ — > oo) behaviour 



4 

X(x,y) ~ (33) 



/(*) ~ ZSTIfMa'M). (32) 
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Thus the iR behaviour of the fermion propagator for an arbitrary value of gauge parameter 
a in our non-local gauge is given by 

s a (^y) t"X-»o> ■ ^ 

\x — y\ t 2 jv 

A special case of eq. (|3*3|) is that the power of|x — y|is3 — 8/ (37r 2 iV) to the leading order 
in 1/iV in the Landau gauge, which is obtained by using a = and «o = 2/3 This 
value was obtained earlier by Aitchison et al j^]. However these authors use a different 
non-local gauge function (the small momentum limit of our gauge function) in the LKF 
transformation equation, and so need to regularize an ultraviolet infinity and also put an 
ultraviolet cutoff scale. Our method is free from these complications. 

Eq. ({B1J) suggests that the iR behaviour for other values of a is again a CFT, albeit 
with a non-zero anomalous dimension for the fermion. We may check this by obtaining the 
dependence of the four-fermion Green function 



(xi,x 2 ;y 1 ,y 2 ) 



8r] lx (x 1 )8t] 72 (x2)5t]s 1 (yi)6r]s 2 (y 2 ) 
on a. From eq. (fTTJj) we get 



(35) 



j=r)=r)=0 



8S 71j y 2 - j s 1 s 2 {xi, X2, yi, y 2 ) [8 (F X1X2 + Fy x y 2 F xi y 2 F X2yi 

—Fxiyi - F X2V2 )\ S 11 , 12; s 1 6 2 (xi,X2;yi,y2)- (36) 

The solution to this equation can be cast in the form 



S^ 2 . Sl52 (x 1 ,x 2 ;yi,y2)- (37) 



' A(gi,g 2 )A(yi,y 2 ) 1 a ~ a ° 

_A(xi, yi)\(x 1 ,y 2 )\(x2, j/i)A(x 2 , 2/2). 

Using eq. ()33|). we may write 

S^, r2 ;6 1 8 2 (xuX2;yi,y2) 

~) S"°^ 2 . SiS2 (x 1 ,x 2 ;yi,y2) (3£ 
for the iR behaviour. Here p and r\ are the conformal invariant cross-ratios 



w 


- x 2 


\Vi 


— 2/2I 


\xi 


-yi\ 


\X2 


— 2/2 | 



\xi — 2/2 1 


\X 2 


— 2/1 1 


W - 2/1 1 


\X 2 


— 2/2 1 



P = i n 1 > V = 1 n r • ( 39 ) 



7 



As S a ° has a structure required by conformal invariance (in the infrared), eq. ()38|) implies 
that S a also has such a structure j7j. This is consistent with an anomalous dimension 
4(a — ao)/(7T 2 N) for the fermion. It is interesting that the angular dependence of the 
scattering amplitude is modified by a simple factor given by a power of p/rj when one 
changes the gauge parameter. 

We now address the gauge dependence of the three-point fermion-photon Green function 



V m ,s{z; x,y) 



(40) 



5j^(z)5r] J (x)Sr] 5 (y) 
From eq. and eq. (|4()jl it follows that 

8V K7tS (z;x,y) = -5F{x,y)V^^(z;x,y) 

+{i/e)d» (5F(z, x) - 8F{z, y)) S jS {x, y). (41) 

For the part of the 3-point function corresponding to a longitudinal photon, we get a simpler 
equation by applying d% on eq. (jJTJ): 

5d^ s (z;x,y) = -6F(x,y)d»V^ >s (z;x,y) (42) 

-ie {6 (d 2 g)-l - 6 (d 2 g)^) S lS (x,y). 

This is satisfied by the Ward identity 

d»V^ s (z;x,y) = -ie [(d 2 g)^ - (d 2 g)^] S 7 s(x,y) (43) 

following from eq. (JSJ). On the other hand for the part V of V that relates to a transverse 
photon, 

/ d z d z v \ 

V^siz; x, y) = ( rj^ - J v »;i,s( z '> x > 3/)> ( 44 ) 

eq. (jlH) leads to 

SV K7! s(z; x, y) = -SF(x, y)V m>s (z; x, y). (45) 

Its solution is 

X(x~y~)) W^^y), (46) 

which is consistent with an anomalous dimension as given in eqs. (I28J) and (J34|) for the 
fermion and a gauge-invariant anomalous dimension for the photon. 

We have shown that the iR behaviour in a class of non-local gauges parametrized by a 
parameter a is given by a CFT with the fermion anomalous dimension depending on the 
parameter a. Using this we now obtain the iR behaviour in the usual class of local gauges. 
This turns out to be very instructive regarding attempts to resum iR divergences of the 
perturbation theory. 

We obtain the fermion propagator in the local gauge corresponding to the Lorentz gauge 
term — (l/(2a)) (OA) 2 by comparing with that for the non-local gauge 

- -L (OA) ■ (l + -j^=] ■ (OA) (47) 



2a v ' V v 73 ^ 
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with the same gauge parameter a. For the local Lorentz gauge 

ae 2 f°° dk ( %m.{k\x - y\) 



2ir 2 J k 2 \ k\x — y\ 

ae 2 , [°° 1 / sinw , 

x - y\ / du— 1 . (4 



2tt 2 J u 2 \ u 

The integral can be evaluated by rewriting the integrand as (u — sinu)(l/u 3 ), integrating 
by parts twice, and using J °° du(smu/u) = ir/2. We thus find F L (x, y) = ae 2 \x — y\/(87i). 
Using eq. (J24*j) the iR behaviour in the local gauge then comes out as 



S L (x,y) ~ exp 



ae 2 



in 



\x - y\ 



*~ 1 (49) 



\x-y 



3 4 "o 



[It may be noted that starting from our non-local gauge, we reach the local Lorentz gauge 
through infinitesimal changes of the gauge function by (formally) varying \x from Ne 2 /8 to 
zero, at a fixed value of a. As fi is decreased, we pass from the \xk 3> k 2 regime to the 
k 2 3> \ik regime, and finally reach F L (x,y) smoothly] 

Now a > for the contribution of the gauge-fixing term to be of the correct sign to 
make the Euclidean functional integral converge. Thus eq. (J49j) tells us that the fermion 
propagator falls off exponentially as if the fermion has developed a mass ae 2 / {8ir)\ However, 
this apparent mass is spurious, since fermion mass cannot be dynamically generated in 
perturbation theory (indeed, the propagator of eq. (J4*9|) is proportional to j> in momentum 
space) . 

This strikingly illustrates the pitfalls in resumming iR divergences in perturbation the- 
ory. There are severe iR divergences in the local Lorentz gauge, because the longitudinal 
part of the photon propagator has a 1/k 2 behaviour in the iR. The cumulative effect is an 
apparent mass term in gauge non-invariant Green functions. The apparent mass is a gauge 
artifact; it does not appear in gauge-invariant correlation functions. 

In this paper, we determined how the iR logarithms of massless QED 3 add up for 
arbitrary values of the gauge parameter in a non-local gauge and also in the usual Lorentz 
gauge. We demonstrated by studying various correlation functions that the iR behaviour 
in the non-local gauge is that of a CFT with a continuously varying anomalous dimension 
for the fermion. We also demonstrated the pitfalls in summing the severe iR divergences 
of the usual Lorentz gauge (the fermion propagator falls off exponentially as if there is a 
fermion mass, which is actually a gauge artifact); thus it is the non-local gauge which is 
suitable for studying the iR behaviour of this theory. The implications of the calculation in 
the non-local gauge for the important issue of anomalous dimension of the gauge-invariant 
dressed fermion will be presented elsewhere jS]. 
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